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Abstract. 

For n- vertex, d-dimensional lattices A with d > 2, the number of spanning trees 
A^st(A) grows asymptotically as exp(7T.ZA) in the thermodynamic limit. We present an 
exact closed-form result for the asymptotic growth constant Zf,cc(d) for spanning trees on 
the d-dimensional body-centered cubic lattice. We also give an exact integral expression 
for Zfcc on the face-centered cubic lattice and an exact closed- form expression for Z488 
on the 4 • 8 • 8 lattice. 
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1. Introduction 

Let G = (V, E) denote a connected graph (without loops) with vertex (site) and edge 
(bond) sets V and E. Let n = v{G) = \V\ be the number of vertices and e{G) = \E\ the 
number of edges in G. A spanning subgraph G' is a subgraph of G with v{G') = \V\, 
and a tree is a connected subgraph with no circuits. A spanning tree is a spanning 
subgraph of G that is a tree (and hence e{G') = n — 1). A problem of fundamental 
interest in mathematics and physics is the enumeration of the number of spanning trees 
on the graph G, Nst{G). This number can be calculated in several ways, including as a 
determinant of the Laplacian matrix of G and as a special case of the Tutte polynomial 
of G m 12] • In this paper we shall present an exact closed-form result for the asymptotic 
growth constant for spanning trees on the d-dimensional body-centered cubic lattice, 
denoted bcc{d), with 6cc(3) = bee. We shall also give an exact integral expression for 
the Zfcc describing the face-centered cubic lattice and an exact closed-form expression 
2:488 for the 4 ■ 8 ■ 8 lattice. A previous study on the enumeration of spanning trees and 
the calculation of their asymptotic growth constants was carried out in Ref . jSj . In that 
work, closed-form integrals for these quantities were given, and from the integral for the 
bcc{d) lattice, an infinite series representation was derived. Our present result for the 
bcc{d) lattice is obtained by summing exactly this infinite series. Similarly, our present 
result for the 4 • 8 ■ 8 lattice is obtained by an exact closed-form evaluation of the integral 
given for this lattice in Ref. 0. 

2. Background and Method 

We briefly recall some definitions and background on spanning trees and the 
calculational method that we use. For G = G(y,E), the degree ki of a vertex Vi & V 
is the number of edges attached to it. A fc-regular graph is a graph with the property 
that each of its vertices has the same degree k. Two vertices are adjacent if they are 
connected by an edge. The adjacency matrix A{G) of G is the nxn matrix with elements 
Aij = 1 if f j and Vj are adjacent and zero otherwise. The Laplacian matrix Q = Q{G) 
is the nxn matrix Q with Qij = ki6ij — Aij. One of the eigenvalues of Q{G) is always 
zero; let us denote the rest as \i{G), 1 < i < n — 1. A basic theorem is that [H |2] 
Nst{G) = (l/n) YYi=i ^i{G). Here we shall focus on fc-regular d-dimensional lattices A. 
For these lattices, if > 2, then in the thermodynamic limit, Nst grows exponentially 
with 77, as n 00; that is, there exists a constant such that Nst{A) ~ exp(?7,ZA) as 
n — i> 00. The constant describing this exponential growth is thus given by 

za = lim n~'^ In Nst{A) . (1) 

where A, when used as a subscript in this manner, implicitly refers to the thermodynamic 
limit of the lattice A. A regular (i-dimensional lattice is comprised of repeated unit cells, 
each containing u vertices. Define a{n, h') as the v x u matrix describing the adjacency 
of the (d-dimensional) vertices of the unit cells n and n', the elements of which are 
given by a{n, n')ij = 1 if t>j G fi is adjacent to Vj G h' and otherwise. Assuming that 
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a given lattice has periodic boundary conditions, and using the resultant translational 
symmetry, we have a(n, h') = a{n — n'), and we can therefore write a{n) = a{hi, ■ ■ ■ , n,])- 
In Ref. jH] a method was derived to calculate Nst{-^) and za in terms of a matrix M 
which is determined by these a{n,n'). For a (i- dimensional lattice, define 

M(^i,---,^,) = A;-l-5^a(n)e*"'^ (2) 

h 

where in this equation 1 is the unit matrix and 6 stands for the (i-dimensional vector 
{ei,---,Od). Then 



/ [UY-]Hdet{M{e„---,eM (3) 

For a fc-regular graph A, a general upper bound is za < Infc. A stronger upper 
bound for a /c-regular graph A with coordination number /c > 3 can be obtained from 
the bound jH E] 

where 

With eq. (0), this then yields [SI 

< In(Cfc) . (6) 

It is of interest to see how close the exact results are to these upper bounds. For this 
purpose, we define the ratio 

where k is the coordination number of A. 
3. hcc{d) Lattice 

For the bcc{d) lattice a unit cell contains utccid) = 2 vertices located at f i = (0, ■ ■ ■ , 0) 
and V2 = (|, ■ ■ ■ , |). This lattice has coordination number fcf,cc(d) = 2'^. Using eq. (jS}, 
Ref. 13] obtained 

Zbcc{d) = d In 2 + hcc{d) (8) 



zk = - 

V 



where 



hcc{d) 



■^-^ j=i j=i 

£[nt]i"(i-n-«.)- M 
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Expanding the logarithm and carrying out the integration term by term yields the 
infinite series representation jS] 

4ccW--2l,^^^5F^J (10) 

We now sum this series exactly. First, 

l^ (£-l)![(2f)!]-^ 

(A;!)2[(2A; + 2)!]'^ 



2^22{^'+i)'^[(A; + l)!]2'i+iA;! 



1 ^ [T{k + l)f[T{2k + ?,)Y 
2d+i 2^ 2(2'=+!)^ [r(A; + 2)]2'i+i A;! ' 



(11) 



Next, we use the duplication formula for the Euler gamma function, 

V{2z) = (27r)-i/2 2^^-^ V{z) V{z + ^) (12) 
with 2; = /c + |, together with r(l/2) = -y/vr, to express 

r(2fc + 3) r(fc + |) 
22'=+ir(fc + 2) r(|) ■ ^ ' 

Substituting this into eq. (fTT|) . we have 

1 - [r(A: + i)p[r(fc + |)/r(|)]-^ 

hcc{d) 2^+1 [r(A; + 2)Y+^ k\ 

= - 2-('^+^) d+2Fd+i([l, 1, 3/2, ■ ■ ■ , 3/2], [2, ■ ■ ■ , 2], 1) (14) 

where there are d + 2 entries in first square bracket [■ ■ ■] and d + 1 entries in the second 
square bracket [■ ■ ■] in the argument, and pFg is the generalized hypergeometric function, 

pF,([ai,---,ap],[6i,---,6g],x) = j (15) 

where c„ = r(c + n)/V{c). Hence, 

= rfln2-2-('='+i) ,+2i^,+i([l,l,3/2,---,3/2],[2,---,2],l) . (16) 

We comment on some special cases. For c? = 1, the 6cc(l) lattice with free 
(periodic) boundary conditions degenerates effectively to a line (circuit) graph, for 
which, respectively, Nst = 1 and Nst = n; in both cases, it follows that Zhcc{i) = 0. 
Using the value 3-F2([l, 1,3/2], [2,2], 1) = 4 In 2, we recover this elementary result. For 
d = 2, the 6cc(2) lattice is equivalent to the square lattice, for which Zgq = (4/7r)/5(2) = 
1.1662436.. pin , where /?(s) = E^=o(-l)"(2^+l)"' and p{2) = C = 0.915965594177.. 
is the Catalan constant. The general result (jH)) with ()14|1 evaluated for d = 2 agrees 
with this, since 4i^3([l, 1, 3/2, 3/2], [2,2,2],1) = 16(ln2 - (2C/7r)). Our general exact 
result for Zbcc{d) provides quite accurate values for higher values of d, which we list in 
Table ^ together with the corresponding ratios (jZj) which give a comparison with the 
upper bound ©. Evidently, the exact values are very close to this upper bound and 
move closer as d increases. 
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Table 1. Values of Zbcc{d) and ncc{d)- 



1 

2 1.166243616123275 0.9587702228064145 

3 1.990191418271941 0.9912457055306051 

4 2.732957535477362 0.9977098978275579 

5 3.447331914522398 0.9993413280070963 

6 4.150116933352462 0.9998002121159708 

7 4.847789269805724 0.9999373061649456 

8 5.543104959793989 0.9999798500846987 

9 6.237305017795394 0.9999934053622532 

10 6.930967870288660 0.9999978103135475 



'2' 2''' V2' '2''' "^^^ V2' 2' 



4. fee Lattice 

The face-centered cubic (fee) lattice has coordination number kfcc 
consisting of the Ufcc = 4 vertices (0,0,0), (0 ^ ^ 
lattice, M {01, 02, 03) is 

/ 12 -(t;2t;3)* -(t^it^s)* -(t^it^2)* \ 

-V2V3 12 -■^^'^2 -'^^'^3 

— ViV^ — flf2 12 ^^2^3 

y —V1V2 —VlV^ ~V2V^ 12 ^ 

where Vj = 1 + e*^^ , j = 1, 2, 3. The evaluation of the determinant yields 

Zfcc = H^2) + - 1^ 1^ —lnF{0i,02,03) 
where, with the abbreviation Cj = cos{0j/2), 



12 and a unit cell 
0). For this 



M(^l,^2,^^3 



(17) 



(18) 



F{01, 02 J 03 



1 xU 1/ 

1 + -^{-C2C3 + C3C1 + CiC2)J 1^1 + -(C2C3 - C3C1 + C1C2J 

1 + ^(C2C3 + C3C1 - C1C2) 1 - ^(C2C3 + C3C1 + C1C2) 



1 - g[(ClC2)' + (C2C3)2 + (c3Ci)2] - — (C1C2C3)' 
^(C1C2C3)^(C? + + C^) + ^[(C1C2)^ + (0203)^ + (030,)'] . 



(19) 

(This corrects an algebraic error in eq. (5.3.3) of Ref. Pj). Evaluating this numerically, 
we find that Zfcc — 2.41292. Substituting Zfcc into eq. ((Tj), we get rj^c — 0.98915, so 
that the upper bound ((Zj) is very close to the actual value of the growth constant. 
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5. 4 ■ 8 • 8 Lattice 

An Archimedean lattice is a uniform tiling of the plane by regular polygons in which 
all vertices are equivalent. Such a lattice can be defined by the ordered sequence of 
polygons that one traverses in making a complete circuit around the local neighborhood 
of any vertex. This is indicated by the notation A = (HiPD' meaning that in this 
circuit, a regular pj-sided polygon occurs contiguously times. We consider here the 
4 ■ 8 • 8 lattice involving the tiling of the plane by squares and octagons. In eq. (4.11) 
of Ref. 1^, the asymptotic growth constant for this lattice was calculated to be 

1 1 d9 dO r 

2:488 = -ln2 + - / / In 7 — 3(cos6'i + cos6'2) — cos6'i cos^2 

2 4 27r 2n L J 

= -ln2 + — / d6\u. 7 — 3 cos ^ + 4sin(^/2)V5-cos0 (20) 

*J 

where the integral on the second line of eq. ()2()|1 is obtained by doing one of the 
two integrations in the expression on the first line. These integrals were evaluated 
numerically to obtain the result z^^^ = 0.786684(1), where the number in parentheses 
indicates the estimated error in the last digit. 

We have derived an exact closed form expression for this integral. We begin by 
recasting the integral in the equivalent form. 



^488 = - In 2 + — / de \n(2 sin(^/2) + Ja + 2sm'^{e/2) 
A 2tx Jq \ V 

ln(^sin(0) + y^l + (l/2)sin2(0) ) . (21) 



That is, 



where 



3 ^ ^7r/2 

-ln2 + - 

4 TT 



Z488 = 7ln2 + 1(1/72) (22) 



4 



1 



TT 



7r/2 







/(«) = -/ rf0 In sin0+ Wl + a^sin^^ ) . (23) 



In eq. (jSHI), with no loss of generality, we take a to be nonnegative. We will give a 
general result for /(a) and then specialize to our case a = l/V^. First, we note that 
/(I) = C/tt, where C is the Catalan constant. Next, assume < a < 1. Taking the 
derivative with respect to a and doing the integral over in eq. (j2SI), we get 
_ -a/2 + (2/7r)tan-^a 

To calculate /(a), we then use /(a) — /(O) = I'{x)dx and observe that 

1 /■'^/^ In 9 9C 

1(0) = - d(P ln(sin(0) + 1) = -— + — (25) 

TT Jo 2 TT 

We also make use of the integrals 

^ dx = --\n{l - a^) (26) 







'l-x2) 2 
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and 



to obtain 



tan 



X 



'l-x^ 



dx = — 



C 



TT 



In 



1 + a 



+ -T 2 
2 



r/ ^ C 1, 

vr 2 



1 



It. 

TT 



1 + a 



if < a < 1, 



(27) 



(28) 



where Ti2(x) is the inverse tangent integral |Sj 

tan-i 



Ti2(x) 



y 



y 



dy 



(29) 



= x[3F2([l,l/2,l/2],[3/2,3/2],-x^)] . 

(Here the arctangent is taken to lie in the range — 7r/2 < tan~^?/ < vr/2.) Evaluating 
our result for J(a) at a = l/-\/2 and substituting into eq. (j221), we obtain the exact, 
closed-form expression 



2^488 



712(3 + 2^2) . 



(30) 



7 + ^Mv^-i) + ^ 

71 Z TT 

The numerical evaluation of eq. (|Hn|) agrees with the evaluation given in Ref. ^ to 
the accuracy quoted there and allows one to obtain higher accuracy; for example, to 
15 significant figures, 2:488 = 0.786684275378832. We note that the Ti2 function also 
appears at intermediate stages in the derivation of Ztri for the triangular lattice . For 
completeness, we have also calculated I{a) for a > 1 with the result 

(a + 1)^1 1„. /I + a' 



TT 2 



2(a- 1) 



It, 

TT 



if a> 1 . 



(31) 
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